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Ultrasonic  consolidation  is  a  solid-state  bonding  process  capable  of  producing  metal  and  metal  matrix 
composite  parts.  In  this  work  a  friction-based  heat  generation  model  is  proposed  to  characterize  the  ther¬ 
mal  development  of  ultrasonically  consolidated  aluminum  foils  and  continuous  fiber  alumina  reinforced 
aluminum  metal  matrix  composite  tape  as  a  function  of  process  control  parameters.  The  friction  coeffi¬ 
cient  between  mating  surfaces  is  determined  experimentally,  and  the  credibility  of  using  both  a  constant 
friction  coefficient  and  a  process  dependent  friction  coefficient  is  assessed.  In  most  cases  a  constant  fric¬ 
tion  coefficient  is  capable  of  producing  results  that  are  within  15%  error;  while  a  process  dependent 
friction  coefficient  achieves  an  average  error  of  7%. 
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1.  Introduction 

Ultrasonic  consolidation  (UC)  is  a  low  temperature  bonding 
process  that  can  be  used  in  the  fabrication  of  metal  and  metal 
matrix  composite  (MMC)  parts.  Layers  of  metallic  films  (foils  or 
prepreg  MMC  tape)  are  built  up  on  top  of  each  other  forming 
the  desired  shape;  this  is  synonymous  to  the  methodology  used 
in  automated  tape  placement  (ATP)  currently  employed  in  the 
thermoplastic  tape  lamination  industry  in  use  by  aerospace  com¬ 
panies.  Tierney  and  Gillespie  (2006)  investigated  and  modeled  the 
in  situ  strength  development  of  the  thermoplastic  ATP  process. 
By  building  up  thin  successive  layers  and  avoiding  bulk  heat¬ 
ing,  residual  stresses  from  fabrication  are  reduced.  UC  provides 
many  advantages  over  traditional  liquid  processing  methods  typ¬ 
ically  employed  in  the  manufacture  of  composites  as  discussed 
by  Doumanidis  and  Gao  (2004).  Weld  temperatures  are  typically 
below  50%  of  the  melt  temperature  and  the  time  at  this  temperature 
is  very  short  (0.02-0.72  s),  thus  allowing  materials  to  retain  much 
of  their  preprocessed  crystallography.  Upon  inspection  of  the  post 
weld  interface  there  is  often  no  indication  of  melting  or  recrystal¬ 
lization  for  low  temperature  welds;  this  has  been  observed  by  both 
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Yang  et  al.  (2009)  and  Clews  (2009).  However,  Clews  does  state  that 
microstructural  changes  could  readily  occur  given  the  right  com¬ 
bination  of  temperature  and  time.  Such  is  likely  the  case  for  welds 
made  by  Mariani  and  Ghassemieh  (2010).  Mariani  and  Ghassemieh 
welded  at  speeds  of  34.5  mm/s  and  predicted  high  temperatures  to 
be  50-80%  of  the  melt  temperature  (330-528  °C).  For  these  hotter 
foil-foil  welds  some  microstructural  changes  were  visible  as  Mari¬ 
ani  and  Ghassemieh  observed  recrystallization  in  a  thin  layer  (5  p,m 
band)  at  the  weld  interface,  between  100  p,m  foils,  using  electron 
backscatter  diffraction.  Largely  retained  crystallography  and  min¬ 
imal  residual  stress  buildups  introduced  by  tape  lamination  may 
reduce  or  eliminate  the  need  for  costly  post-processing  heat  treat¬ 
ments.  Localized  heating,  versus  bulk  heating,  also  helps  reduce 
processing  costs  by  decreasing  the  amount  of  energy  required  dur¬ 
ing  fabrication.  Additionally,  low  temperature  layered  processing 
facilitates  the  use  and  placement  of  thermally  sensitive  materials 
(e.g.  embedded  sensors  and  fibers)  that  could  be  otherwise  dam¬ 
aged  or  imprecisely  placed  using  liquid  processing.  This  has  been 
demonstrated  by  various  researchers.  Cheng  et  al.  (2007)  success¬ 
fully  embedded  nickel  based  thin  film  thermocouples  into  copper 
work  pieces  via  UC.  Siggard  et  al.  (2006)  consolidated  USB-based 
sensors  into  aluminum.  Kong  et  al.  (2004b)  embedded  shape  mem¬ 
ory  alloy  (SMA)  fibers  into  aluminum  to  create  adaptive  structures 
for  aerospace  applications.  After  embedding  SMA  fibers,  Kong  and 
Soar  built  upon  their  previous  work  by  using  SiC  (2005a)  for  struc¬ 
tural  applications  and  optical  fibers  (2005b)  for  data  transport. 
With  UC  it  is  also  possible  to  weld  a  variety  of  dissimilar  mate¬ 
rials  including:  aluminum,  brass,  stainless  steel,  super  alloys,  SiC 
fiber,  and  MMC  tape.  Janaki  Ram  et  al.  (2007)  consolidated  many  of 
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these  material  combinations  and  investigated  the  bonding  through 
microscopy. 

Fig.  1  is  a  schematic  of  the  key  elements  involved  in  UC. 

The  cylindrical  horn  (also  referred  to  as  the  sonotrode )  applies 
a  normal  force,  F,  which  brings  the  top  tape  in  contact  with  the 
material  to  which  it  will  be  bonded.  Sufficient  force  is  required  to 
ensure  intimate  contact  at  the  interface.  Processing  time,  t,  is  con¬ 
trolled  by  the  contact  length  (the  longitudinal  distance  over  which 
force  is  applied),  lc,  and  linear  weld  speed  (tangential  speed  at  the 
horn’s  surface,  determined  via  the  horn’s  RPM  and  diameter),  s.  s 
and  t  are  defined  by  Eqs.  (1 )  and  (2). 


tJj  (2) 

The  horn  oscillates  at  a  fixed  frequency  of  20  kHz,  /,  and  peak- 
to-peak  amplitude,  %.  The  number  of  oscillation  cycles,  N,  a  sample 
is  subjected  to  for  any  given  weld  is  determined  by: 

N  =  tf  =  y  (3) 

The  surface  of  the  horn  is  knurled.  This  provides  a  firm  grip 
between  the  horn  and  the  upper  tape  thus  preventing  slip,  and 
consequently  welding,  at  the  tool-tape  interface.  Thereby  relative 
motion  occurs  between  the  upper  tape  and  the  substrate  result¬ 
ing  in  frictional  work.  Friction  causes  abrasion  of  the  contacting 
surfaces  smoothing  out  irregularities,  breakup  and  dispersion  of 
surface  oxides  and  asperities,  and  heat  generation  from  dissipated 
frictional  work  which  all  promote  bonding  and  welding.  Additional 
details  regarding  the  specifications  of  the  welder  used  in  this  study 
as  well  as  the  materials  processed  will  be  given  in  the  equipment 
and  materials  section. 

The  exact  bonding  mechanisms  involved  in  UC  are  not  well 
defined.  It  has  been  suggested  that  diffusion  and  plastic  deforma¬ 
tion  aid  in  the  bonding  process  as  theorized  by  Neppiras  (1965) 
and  discussed  by  Hazlett  and  Ambekar  (1970).  Janaki  Ram  et  al. 
(2007)  attributed  flow  lines  to  plastic  deformation  which  permitted 
the  embedment  of  fibers  during  UC.  Yang  et  al.  (2009)  has  uti¬ 
lized  various  optical  techniques  (orientation  imaging  microscopy, 
OIM,  and  X-ray  energy  dispersive  spectroscopy,  EDS)  that  quanti¬ 
fied  the  presence  of  plastic  deformation  and  diffusion  during  UC.  It 
is  also  known  that  the  bonding  mechanisms  of  ultrasonic  welding 
are  heat-assisted.  Though  not  required,  the  application  of  addi¬ 
tional  heat  (pre-heating)  facilitates  the  welding  process  as  reported 
by  Neppiras  (1965).  Plastic  deformation  occurs  from  friction,  the 
applied  load,  and  thermal  and  acoustic  softening.  Acoustic  soft¬ 
ening  can  result  when  a  material  is  subjected  to  high  frequency 
loads.  Langenecker  (1966)  showed  that  the  effects  of  ultrasonic 
exposure  has  much  the  same  effect  as  elevated  temperatures  do 
on  the  yield  strength  through  the  comparison  of  stress  vs.  elon¬ 
gation  plots  at  varying  temperature  and  ultrasonic  energy  levels. 
In  turn,  this  makes  localized  plastic  deformation  occur  more  read¬ 
ily  which  is  required  to  bring  mating  surfaces  together  and  fill  the 
gaps.  Gaps  can  be  naturally  occurring  due  to  material  variation  or 
induced  by  embedding  materials  (e.g.  fibers,  sensors,  etc.).  While 
some  material  mixing  may  occur  providing  a  mechanical  bond  from 
plastic  flow,  the  experimentally  measured  concentration  gradients 
obtained  through  EDS  by  Gunduz  et  al.  (2005)  and  Mariani  and 
Ghassemieh  (2010)  suggest  diffusion  occurs  along  the  boundaries 
between  materials  producing  superior  atomic  bonds.  However, 
since  surfaces  must  be  in  close  contact,  the  applied  pressure  and 
plastic  deformation  aid  diffusion.  Diffusion  is  a  temperature  depen¬ 
dent  process,  as  is  yield  strength,  thus  the  processing  temperature 
and  weld  time  are  important  factors  in  the  bonding  process. 


Temperature  increases  are  attributed  to  two  sources:  bulk  plas¬ 
tic  deformation,  and  interfacial  friction.  The  latter  being  the  focus 
of  this  paper.  It  is  assumed  that  slip  occurs  at  the  interface  between 
materials  being  consolidated  for  the  entire  duration  of  the  weld¬ 
ing  process,  which  promotes  friction  and  deformation  of  material 
at  the  weld’s  interface.  For  this  study  thermal  contributions  from 
deformation  are  neglected  and  only  frictional  contributions  are 
accounted  for. 

1.1.  Deformational  heat 

Thermal  contributions  from  deformation  are  negligible  for  sev¬ 
eral  reasons;  though,  this  is  not  immediately  obvious.  To  assess  the 
amount  of  heat  generated  by  deformation  the  amount  of  plastic 
work  needs  to  be  determined.  Plastic  work  per  unit  volume  can 
be  described  by  the  product  of  the  yield  stress  and  the  change 
in  plastic  strain  developed.  If  the  entire  sample  flowed  plastically 
without  slipping,  and  all  plastic  work  was  dissipated  as  heat,  defor¬ 
mational  heat  generation  could  produce  approximately  five  times 
the  energy  of  the  frictional  contribution  presented  in  this  paper  (see 
Section  2.1  for  heat  generation  flux  derivation),  assuming  normal 
room  temperature  properties.  However,  numerous  factors  limit 
and  diminish  the  ability  for  deformation  to  contribute  significantly 
to  the  thermal  development  during  UC  of  foils  and  tapes.  Small 
temperature  gradients  through  the  material’s  thickness  will  cause 
a  gradient  of  the  yield  stress.  The  weld’s  interface  temperature  will 
be  the  highest  (either  from  localized  plastic  deformation  and/or 
friction).  Areas  of  highest  temperature  will  have  the  lowest  yield 
strength.  By  accounting  for  a  small  amount  of  thermal  variation  it 
is  expected  that  only  material  near  the  interface  will  yield  plas¬ 
tically.  This  has  been  confirmed  experimentally.  The  entire  bulk 
material  does  not  flow  plastically.  In  a  technical  report  promot¬ 
ing  UC  White  (2003)  stated  that  Plastic  flow  is  confined  to  a  thin 
interfacial  layer  10-20  |xm  thick,  which  is  at  most  10%  of  the  sam¬ 
ple’s  thickness.  This  was  proven  through  grain  boundary  mapping 
of  the  cross-section  of  a  consolidated  sample  by  means  of  electron 
backscatter  diffraction  analysis  by  Mariani  and  Ghassemieh  (2010) 
in  their  microstructural  investigation  of  UC  Al.  For  the  case  of  MMC 
tape  in  this  study  only  the  matrix  will  accumulate  plastic  strain 
as  the  brittle  ceramic  fiber  reinforcement  would  break  if  yielded 
plastically.  Thus  the  volume  of  material  behaving  plastically  is  fur¬ 
ther  reduced  via  the  fiber  volume  fraction.  Inside  the  plastic  zone, 
only  33%  of  the  plastic  work  is  expected  to  be  converted  into  heat 
for  UC  aluminum  (for  non-UC  plastic  deformation  typically  90% 
conversion  is  assumed).  The  remaining  energy  is  stored  in  the  mate¬ 
rial’s  microstructure.  This  assumption  is  also  used  in  Zhang  and 
Li’s  (2009)  dynamic  thermal-mechanical  model  of  the  first  0.0025s 
of  an  ultrasonic  weld  and  is  based  upon  the  work  by  Hodowany 
(1997)  and  Ravichandran  et  al.  (2002).  Hodowany  (1997)  and 
Ravichandran  et  al.  (2002)  carried  out  experiments  measuring  the 
percentage  of  work  dissipated  as  heat  as  a  function  of  the  plastic 
strain  developed.  As  a  consequence  of  the  aforementioned  factors 
present  in  UC  the  potential  contribution  of  deformation  is  reduced 
from  458%  of  the  frictional  contribution  to  6%.  Further  reduction  of 
the  deformational  heat  contribution  results  from  the  decrease  in 
material  strength.  As  stated  in  Hibbeler’s  (2002)  textbook,  using  a 
perfectly  plastic  material  model,  the  stress  level  performing  work 
will  be  limited  by  the  materials  yield  strength.  The  experimen¬ 
tal  results  shown  in  Langenecker’s  (1966)  tensile  testing  indicated 
both  elevated  temperatures  and  the  presence  of  ultrasonic  energy 
significantly  decreases  the  yield  strength  of  Aluminum.  For  UC  con¬ 
ditions  similar  to  this  study  the  reduction  in  yield  stress  from  mild 
heating  (100°C)  and  typical  ultrasonic  softening  (20,000W/m2)  is 
on  the  order  of  80%  as  indicated  by  Siddiq  and  Ghassemieh  (2009) 
in  their  theoretical  analysis  investigating  the  thermal  and  acous¬ 
tic  softening  exhibited  during  UC.  Consequently,  deformational 
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Fig.  1.  UC  Schematic:  left,  isometric  view  of  a  tape  being  welded  to  a  substrate;  right,  through-thickness  cross-section  of  tape  and  substrate. 


contributions  are  reduced  to  approximately  1%  (for  MMC  tape) 
of  the  frictional  contribution  using  this  simplified  calculation.  For 
6061 -T6  foils  this  increases  to  7%,  and  for  3003-H18,  4%.  In  2009, 
Zhang  and  Li  proposed  a  more  detailed  theoretical  dynamic  FEA 
simulation  of  the  UC  process  which  also  supports  the  conclusion 
that  deformational  energy  is  negligible;  after  the  20th  oscillatory 
cycle  (0.001  s)  the  plastic  heat  generation  rate  saturated  and  con¬ 
tributed  to  only  0.3%  (for  3003-H1 8,  ay  =  1 86  MPa)  of  any  additional 
heat  accumulated  at  the  weld’s  interface.  The  back  of  the  enve¬ 
lope  calculation  performed  in  this  section  assumes  uniform  stress 
across  the  width  (vs.  three  dimensional)  and  does  not  allow  slip; 
consequently,  it  is  a  very  conservative,  upper  limit,  approach  to 
approximating  the  potential  contribution  of  deformational  heat 
relative  to  frictional  heating. 

Since  deformation  is  a  negligible  thermal  contributor  the  objec¬ 
tive  of  this  work  is  to  explore  the  role  of  friction  in  increasing  the 
interface  temperature  and  its  dependence  on  the  machine  variables 
such  as  speed,  amplitude  and  applied  force.  Section  2  will  pro¬ 
pose  a  model  to  relate  temperature  through  energy  balance  to  the 
friction  coefficient.  The  value  of  the  friction  coefficient  will  be  deter¬ 
mined  by  conducting  experiments  under  various  speeds,  applied 
forces  and  amplitudes  on  an  ultrasonic  welder  and  by  matching 
the  experimental  temperatures  measured  with  an  1R  camera  with 
model  predictions.  Some  of  the  challenges  in  this  characterization 
technique  will  be  outlined. 

f  .2.  Friction 

Friction  is  the  resulting  force  that  opposes  motion  between  two 
contacting  surfaces.  In  Bowden  and  Tabor’s  (1950)  book,  much 
work,  experimentally  and  theoretically,  was  done  to  advance  the 
understanding  of  the  mechanisms  involved  in  the  friction  and 
lubrication  of  solids.  They  report  that  the  classical  expression  quan¬ 
tifying  friction,  Eq.  (4),  dates  back  as  far  as  Leonardo  da  Vinci,  but 
most  credit  it  to  the  work  of  Amontons  and  Coulomb.  The  law 
states  that  the  force  of  sliding  friction,  Ffr,  is  proportional  to  its  nor¬ 
mal  force,  Fn.  The  constant  of  proportionality,  /x,  is  interchangeably 
referred  to  as  the  coefficient  of  friction  or  the  friction  coefficient. 

Ffr  =  ri-FN  (4) 

J.2.  J.  Friction  coefficient  variability  and  trends 

The  friction  coefficient  is  a  dimensionless  parameter  deter¬ 
mined  experimentally  that  quantifies  the  effects  of  a  variety  of 
elements  that  influence  the  resistance  to  slip  between  materi¬ 
als.  Blau  (2009)  wrote  a  more  modern  book  referencing  much  of 
the  more  modern  studies  and  standards  pertaining  to  friction.  An 
underlying  theme  of  the  book  is  that  beyond  material  properties 


and  surface  characteristics  there  is  a  plethora  of  factors  that  influ¬ 
ence  the  coefficient  of  friction,  several  of  which  vary  during  UC. 
The  input  parameters  involved  in  welding  are  clamping  force,  F, 
oscillation  amplitude,  X,  and  weld  speed,  s.  All  of  these,  as  well  as 
temperature,  have  been  shown  to  influence  the  friction  coefficient. 

Clamping  force,  F,  is  proportionally  related  to  the  applied  pres¬ 
sure  through  the  contact  area,  lcw,  so  any  influence  pressure  may 
have  on  the  friction  coefficient  will  be  evinced  by  variations  in  F. 
In  cyclic  fretting  fatigue  testing  done  by  Naidu  and  Raman  (2005), 
applied  slip  displacement  is  prescribed  via  the  applied  cyclic  stress. 
Applied  stress  is  proportional  to  strain  (cr  =  £e),  and  strain  is  pro¬ 
portional  to  the  change  in  gauge  length  (e  =  A L/L).  When  a  fretting 
pad  is  applied  to  the  elongating  bar,  interfacial  slip,  X,  will  also 
be  proportional  to  AL.  Therefore,  the  induced  slip  displacement  is 
proportional  to  applied  cyclic  stress.  Thus  variations  in  the  friction 
coefficient  observed  by  Naidu  and  Raman  (2005)  as  a  function  of 
cyclic  stress  can  be  similarly  related  to  changes  in  X.  This  approach 
is  taken  by  Siddiq  and  Ghassemieh’s  (2008)  in  their  UC  study  to 
determine  /x  as  a  function  of  weld  amplitude,  X.  The  third  input 
parameter,  weld  speed,  s,  is  inversely  proportional  to  N,  the  num¬ 
ber  of  cycles  as  shown  in  Eq.  (3).  In  summary,  F  is  proportional  to 
contact  pressure,  X  is  proportional  to  cyclic  stress,  and  s  is  inversely 
proportional  to  N. 

The  empirical  trends  of  Naidu  and  Raman's  testing  as  related  to 
the  effect  of  contact  pressure,  cyclic  stress,  and  N  on  /j,  is  presented 
in  Fig.  2  for  reference.  In  2007  Zhang  and  Li  submitted  a  conference 
paper  detailing  an  earlier  version  of  their  numerical  model  of  the 
UC  process.  Part  of  this  work  included  investigating  the  influence 
of  temperature  on  the  sliding  friction  coefficient  of  self  mated  Al. 
In  Fig.  2  the  data  from  Zhang  and  Li’s  (2007)  study  on  /x  vs.  T  is 
presented  for  reference. 

Through  inference  the  influence  of  X,  F,  and  s  is  not  insignificant. 
As  a  result,  a  single  valued  constant  friction  coefficient,  /x constant , 
may  not  be  sufficiently  accurate  in  predicting  the  heat  generated 
during  the  process.  Therefore,  a  coefficient  of  friction,  /J-RSM’  which 
depends  on  X,  F,  and  s,  will  also  be  investigated.  This  work  will 
explore  this  dependence  and  compare  the  error  introduced  due  to 
use  of  a  constant  friction  coefficient. 


2.  Thermal  friction  model 

Two  different  models  for  predicting  the  processing  temperature 
during  UC  will  be  employed  in  this  study.  Both  utilize  the  same  heat 
flux  term  based  on  frictional  work.  However,  the  first  model  is  ana¬ 
lytic  and  spatially  invariant,  the  other  is  a  numerical  finite  element 
analysis  (FEA)  two-dimensional  model.  While  the  FEA  solution  is 
more  detailed  and  can  account  for  spatial  variations,  it  requires 
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Fig.  2.  Top  and  lower  left,  Naidu  and  Raman’s  (2005)  friction  trends  as  determined  through  uniaxial  cyclic  fretting  fatigue  testing.  Lower  right,  Zhang  and  Li  (2007)  experimental 
ji  vs.  T  results  for  Al-Al  sliding  contact  (plots  reproduced  with  permission]. 


more  time  to  obtain  results  than  the  analytical  solution.  Therefore, 
if  sample  properties  allow  (e.g.  low  Biot’s  number;  high  internal 
conduction,  and  low  convective  loss,  h,  and  small  sample  size,  1), 
the  analytic  approach  may  be  more  appealing  in  an  environment 
where  FEA  software  is  not  easily  utilized.  This  is  further  explained 
in  Section  2.2. 


2.2.  Frictional  heat  generation  term 

Friction  is  a  non-conservative  force.  Assuming  all  work  done  by 
friction  is  dissipated  as  heat  over  the  contact  area  forms  the  basis  of 
the  friction  model.  The  rate  at  which  work  is  done  (i.e.  power)  is  the 
product  of  the  force  of  friction  and  the  average  oscillatory  speed,  v. 
Dividing  the  friction  power  by  the  contact  area,  w  ■  lc,  defines  the 
frictional  heat  flux.  This  frictional  flux  can  be  written  as: 


% 


Ffr v 
wlc 


(5) 


The  average  oscillatory  speed  is  dependent  on  frequency,/,  and 
amplitude.  Modeling  the  oscillatory  motion,  A,  as  a  sinusoidal  wave 
results  in  Eq.  (6). 


sin  2n ft 


(6) 


Differentiating  A  with  respect  to  time,  f,  yields  the  velocity.  The 
average  oscillatory  speed  is  the  average  magnitude  of  the  veloc¬ 
ity.  This  is  determined  by  integrating  the  absolute  value  of  velocity 
over  one  period  and  multiplying  it  by  the  frequency.  The  resulting 
expression  is: 


v  = 


AA 

dt 


2Xf 


(7) 


By  combining  Eqs.  (4),  (5)  and  (7)  with  F^  =  F,  the  horn  applied 
clamping  force,  the  frictional  heat  generation  flux  can  be  written 
as: 


n„  ri F2Xf 

wlc 


(8) 


2.2.  Model  implementation 

The  energy  balance  applied  to  the  region  of  interest  will  allow 
one  to  solve  for  the  temperature  field  as  a  function  of  time  and  posi¬ 
tion.  The  friction  is  generated  at  the  interface  so  it  can  be  introduced 
as  a  flux  boundary  condition.  Before  we  embark  on  solving  the 
three  dimensional  time  dependent  problem,  it  may  be  worthwhile 
to  check  if  one  can  simplify  the  problem  to  just  time  dependent 
by  evaluating  the  Biot  number.  When  the  Biot  number  is  much 
less  than  one,  there  is  usually  not  much  of  temperature  gradient  in 
the  domain  of  interest  and  one  can  use  lumped  parameter  analysis 
and  solve  analytically  for  the  temperature  as  a  function  of  time.  In 
Incropera  and  DeWitt’s  (2001 )  introductory  heat  transfer  textbook 
the  Biot  number  is  defined  as: 


The  Biot  number  is  the  ratio  of  convective  losses  to  internal 
conduction.  When  materials  are  thin  (small  l )  or  highly  conduc¬ 
tive  (large  k),  the  internal  temperature  variation  is  small.  Thus,  for 
Bi<0.1  (thermally  thin)  Incropera  and  DeWitt  (2001)  state  that 
a  material’s  internal  temperature  can  usually  be  assumed  to  be 
spatially  invariant  making  analytic  temperature  solutions  possible. 
Conversely,  thick  samples  with  low  conductivity  and  large  con¬ 
vective  losses  will  have  spatially  variant  temperatures  (thermally 
thick),  which  in  turn  are  more  involved  to  solve  analytically.  The 
Biot  number  is  about  0.04  for  welded  foils  and  0.08  for  an  MMC 
tape  welded  to  an  Al  substrate.  Both  tape  and  foils  have  low  Biot’s 
numbers,  but  since  foils  and  tapes  have  high  aspect  numbers  (i.e. 
the  ratio  of  cross-sectional  length  to  thickness),  the  amount  of 
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Fig.  3.  Temperature  variation  for  a  typical  weld:  (a)  through  the  tape’s/foils'  thickness  at  the  midplane,  (b)  across  the  width  (from  the  sample’s  edge  to  the  midplane)  at  the 
weld  interface.  Note  scale  differences  for  temperature  and  position  for  each  plot. 


thermal  variation  in  both  the  thru  thickness  direction  and  across 
the  width  at  the  weld  interface  should  be  investigated.  Using  the 
finite  element  analysis  approach  detailed  later,  the  following  plots 
were  generated  illustrating  the  amount  of  temperature  variation 
for  a  typical  weld  through  the  tape’s/foils’  thickness  at  the  mid¬ 
plane,  Fig.  3(a),  and  across  the  width  (from  the  sample’s  edge  to  the 
midplane)  at  the  weld  interface,  Fig.  3(b).  The  average  interfacial 
temperature  for  the  tape  and  foils  is  the  same,  102  °C. 

In  the  thickness  direction  both  foils  and  tapes  show  little  varia¬ 
tion  in  temperature.  The  coefficient  of  variation  (CV)  is  1%  for  foils 
and  MMC  tape.  Flowever,  across  the  width  larger  variations  are 
observed.  The  CV  across  the  width  is  3%  for  foils  and  10%  for  tape. 
Therefore,  the  foils  used  in  this  study  will  be  considered  thermally 
thin,  whereas  the  MMC  tapes  are  thermally  thick.  Therefore,  the 
solution  approach  will  differ  depending  on  the  materials  being  con¬ 
solidated.  Foil  temperatures  will  be  solved  using  lumped  parameter 
analysis,  while  tape  temperatures  will  be  predicted  with  a  commer¬ 
cial  FEA  software  package,  COMSOL  Multiphysics  (2008). 

2.2.1.  Lumped  parameter  analysis 

Since  the  temperature  of  welded  foils  does  not  vary  spatially, 
temperature  in  the  foil  changes  only  with  time.  The  first  step  is 
to  perform  an  energy  balance  over  the  control  volume  of  interest, 
the  foil  contact  area  and  thickness.  Any  energy  put  into  the  control 
volume  that  is  not  dissipated  into  the  surroundings  is  stored  as 
described  in  Eq.  (10). 

Ejn  —  Eout  =  ^stored  =  PCpV  (10) 

Fig.  4  illustrates  the  control  volume  for  two  identical  foils  being 
welded  together.  Thermal  losses  into  the  air,  horn,  and  supporting 
anvil  (flat  knurled  surface  that  secures  the  lower  foil  preventing 
slip  during  UC)  are  noted  with  outward  arrows.  The  center  shaded 
area  represents  the  friction  heat  flux,  which  is  the  heat  input 
into  the  control  volume. 


Fig.  4.  Foil-foil  control  volume  and  thermal  boundary  conditions. 


Using  Eq.  (10)  in  conjunction  with  the  thermal  loads  from  Fig.  4 
the  following  equation  is  obtained: 

f)Cp V =  Ejn  —  E0ut  —  Qjr  ■  tv/c  —  (hroom^dlc  +  djlornwlc  -F  hanvjjWlc ) 
x(T  —  Too)  (11) 


By  assuming  pt  and  the  convection  coefficients  are  independent 
of  temperature  and  time,  Eq.  (11),  which  is  a  first  order  ODE,  can 
be  integrated  to  yield 


T  Too+  H(\-  e(-Ht/f>cp'/)) 


where  the  volume,  17;  time,  t;  and  equivalent  convective  losses,  H ; 
are  given  by: 

V  =  2  dwlc 

(13) 

t  =  k 

S 

(14) 

H  =  (hroom^dlc  +  hhomWlc  +  hanvnWlc) 

(15) 

While  losses  to  air,  firoom,  can  be  prescribed  based  on  typical 
metal-to-still  air  values,  and  hanvu  cannot  be  predetermined. 
Thus,  H  is  found  empirically  through  inverse  modeling  using  a  two- 
dimensional  FEA  model  and  Eq.  (12).  This  allows  conductive  losses 
to  be  effectively  converted  into  a  convective  loss  term.  The  result¬ 
ing  values  for  H  are  between  3.2  and  3.9  W/K  depending  on  the 
total  weld  time.  Since  convective  coefficients  were  assumed  to  be 
independent  of  time  during  integration,  solutions  to  Eq.  (12)  will 
be  slightly  off.  For  the  purpose  of  this  study  this  is  inconsequen¬ 
tial  since  the  temperature  at  the  end  of  the  weld  is  what  will  be 
predicted  and  measured. 

2.2.2.  Finite  element  analysis 

The  FEA  analysis  does  not  provide  the  level  of  insight  of  an 
analytic  solution,  but  fewer  assumptions  are  required,  compli¬ 
cated  geometries  can  be  modeled,  and  graphical  representation  of 
results  can  be  easily  obtained.  It  is  applied  to  the  MMC  tape  as 
the  Biot’s  number  is  close  to  0.1  and  it  was  shown  that  there  is 
spatial  variation  in  the  temperature  field  (Fig.  3(b)).  Fig.  5  shows 
the  through-thickness  two-dimensional  FEA  model:  the  bound¬ 
ary  conditions,  mesh  properties,  and  heat  flux  boundary  condition. 
Due  to  the  confined  space  between  the  horn  and  substrate  a  low 
convective  loss  was  assumed,  h  =  5  W/m2  K,  and  applied  on  the 
dashed  boundaries.  The  initial  temperature  is  taken  as  ambient 
temperature,  Too.  A  reduced  half-symmetry  geometry  was  used  to 
minimize  computational  resources,  while  still  being  large  enough 
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that  thermal  gradients  are  zero  on  the  fixed  temperature  bound¬ 
aries  (indicated  by  a  bold  line).  Mesh  refinement  was  performed  to 
assure  convergence.  For  the  tape-substrate  FEA  model  a  single  tape 
replaces  the  foils,  a  substrate  replaces  the  anvil,  and  is  applied 
at  the  tape-substrate  interface.  We  can  solve  for  temperature  as  a 
function  of  time  in  the  entire  domain. 

3.  Equipment,  materials  and  experiments 

In  this  section  relevant  laboratory  equipment  will  be  discussed. 
The  geometric  and  thermal  properties  of  the  materials  used  will 
also  be  specified.  All  sample  preparation  and  testing  methods  are 
detailed.  The  parameter  selection  criteria  and  parameter  arrays  are 
defined.  And  the  process  for  obtaining  the  friction  coefficient  is 
explained. 

3.1.  Equipment 

The  primary  equipment  used  in  the  study  was  an  infrared  cam¬ 
era  and  a  customized  ultrasonic  welder.  The  infrared  camera  is 
used  to  measure  temperature  during  the  consolidation  process.  The 
ultrasonic  welder  is  used  to  consolidate  an  MMC  tape  or  metallic 
foil  to  a  metal  substrate  or  another  foil. 

3.  i .  t .  Infrared  camera 

Heat  given  off  by  objects  can  be  seen  in  the  infrared  (IR)  spec¬ 
trum.  To  measure  the  processing  temperature  of  foils  and  tapes  a 
FL1R  Thermovision  Alert,  model  194,  IR  camera  was  used.  The  cam¬ 
era  was  mounted  in  front  of  the  welder  at  a  fixed  distance/angle, 


and  focused  on  the  nip  point  between  the  horn  and  tape  (or 
foil).  Temperature  dependent  emissivities  were  calibrated  using 
the  same  surrounding  geometries  and  camera  positioning  on  a 
hotplate.  Fig.  6  shows  the  welding  direction,  overlaid  IR  image,  tem¬ 
perature  contours  acquired  (across  the  tape  width),  and  camera 
positioning.  The  camera  recorded  in  real  time  during  the  welding 
process.  Each  temperature  contour  captured  corresponds  to  the 
temperature  profile  of  the  nip  point  from  a  single  camera  frame. 
This  temperature  profile  is  that  of  the  sample  just  as  welding  is 
completed. 

3.1.2.  Ultrasonic  welder 

The  welds  were  made  with  a  modified  seam  welder  purchased 
from  AmTech.  The  knurled  Ti-6A1-4V  horn  has  a  diameter  of 
146  mm,  and  can  rotate  at  speeds  up  to  150rpm  resulting  in  weld 
speeds  up  to  1200  mm/s,  via  Eq.  (1).  The  prescribed  amplitude  of 
oscillation  and  fixed  frequency  are  regulated  using  feedback  con¬ 
trol  built  into  the  welder.  The  frequency  used,  20  kHz,  is,  as  reported 
by  Neppiras’  (1965)  investigation  of  the  physical  mechanisms 
involved  in  ultrasonic  welding,  typical  to  ultrasonic  welding.  The 
clamping  force  is  controlled  with  a  pressure  regulated  pneumatic 
cylinder.  The  available  parameter  ranges  (Table  2)  for  the  welder 
are  quite  broad.  This  allows  for  greater  flexibility  and  understand¬ 
ing  of  current  and  future  welding  geometries  and  material  pairings. 
Depending  on  materials  consolidated  (foil-foil  vs.  tape-substrate) 
contact  lengths,  lc,  varied  (see  Table  1).  From  Askeland  and  Phule’s 
(2003)  material  science  and  engineering  textbook  it  can  be  gath¬ 
ered  that  this  is  because  the  matrix  metal  in  the  MMC  tape  is  pure 
Al,  which  has  a  yield  strength  8-16  times  lower  than  the  foils, 


Table  1 

Material  properties. 


Material 

lc  (mm) 

w  (mm) 

d  (mm) 

Cp  O/kg  I<) 

l<  (W/mK) 

P  (kg/m3) 

Foil 

A16061-T6 

2.55 

12.7 

0.1 

896 

167 

2700 

Substrate 

Al  6061 -T6 

7.22 

25.4 

12.7 

896 

167 

2700 

Tape 

A1/A1203  (57%  FVF) 

7.22 

10 

0.36 

802 

81 

3388 
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Fig.  6.  IR  camera  positioning  and  acquired  data. 


Table  2 

Current  welder  parameter  ranges. 


Available  parameter  range 

Amplitude  (p,m) 

6-45 

Force  (N) 

300-5500 

Speed  (mm/s) 

0-1200 

17-34  MPa  vs.  276  MPa.  A  lower  yield  strength  means  that  a  larger 
amount  of  contact  area  is  needed  to  support  the  same  load.  Addi¬ 
tional  detail  regarding  the  contact  length  and  the  procedure  for 
determining  lc  will  be  discussed  further  in  Section  3.3.1. 

3.2.  Materials 

The  foils  and  substrates  used  in  this  study  were  aluminum 
alloy  6061  heat  treated  to  the  T6  condition.  The  Askeland  and 
Phule  (2003)  textbook  reports  that  the  T6  temper  designation  indi¬ 
cates  the  material  has  been  solution  treated  and  artificially  aged 
to  strengthen  it.  The  foils  in  this  study  were  heat  treated  in  an 
oxygenated  environment  resulting  in  a  very  thick  oxide  layer.  This 
layer  can  be  removed  with  a  low  concentration  nitric  acid  bath 
(3%  for  18  h).  Foils  treated  with  nitric  acid  will  be  referred  to  as 
cleaned  foils.  Untreated  foils  will  be  referred  to  as  oxidized  foils. 
This  step  was  not  necessary  for  substrate  welding  as  no  thick  oxide 
existed;  implying  that,  unlike  the  foils,  the  heat  treatment  process 
was  performed  in  a  more  inert  environment.  The  MMC  tapes  were 
provided  by  the  Army  Research  Laboratory.  The  tapes,  commer¬ 
cially  known  as  MetPreg,  were  fabricated  by  Touchstone  Research 
Laboratory  (TRL).  MetPreg  consists  of  Nextel  610  alumina  fibers  in 
a  pure  aluminum  matrix.  The  roll  used  in  this  study  was  57%  fiber 
volume  fraction  (FVF).  Micrographs  of  a  cross-section  of  MetPreg 
are  shown  in  Fig.  7.  There  is  a  large  amount  of  variation  in  thickness 
from  the  edge  (0.43  mm)  of  the  specimen  to  the  center  (0.32  mm), 
30%  difference.  After  UC  tape  thickness,  and  variation,  is  reduced 
(-8%),  and  width  increases  (+9%).  For  thermal  modeling  the  tape 
geometry  was  held  constant  and  assumed  rectangular;  an  in  situ 
geometry  was  used  that  consisted  of  an  average  of  pre  and  post- 
processed  measured  tape  dimensions  (Table  1 ).  Foils  on  the  other 
hand  had  no  measureable  variation  in  thickness  or  width,  pre  or 
post  UC.  Unlike  the  pure  Al  matrix  of  the  MMC  tape,  the  high  yield 
strength  of  Al  6061 -T6  is  more  resistant  to  plastic  deformation.  Thus 
for  foils  the  model  geometry  was  also  the  pre  and  post  weld  mea¬ 
sured  thickness  and  width.  The  fact  that  the  MMC  tape  deforms 
during  UC  will  not  significantly  affect  the  thermal  development, 


since  deformational  contributions  are  expected  to  be  on  the  order 
of  1%,  as  indicated  in  Section  1.1.  Geometric  and  thermal  properties 
for  all  consolidated  samples  are  summarized  in  Table  1. 

3.3.  Experiments 

The  primary  goal  of  the  experiments  is  to  record  the  tempera¬ 
ture  of  the  foil  or  tape  during  the  UC  process.  This  temperature  is 
then  compared  to  the  values  predicted  from  the  model  based  on 
the  best-fit  value  for  a  single  valued  constant  friction  coefficient, 
ftconstant » and  for  a  coefficient  of  friction,  jlrsM'  that  is  allowed  to  vary 
with  the  machine  variables.  Once  ft  is  known,  a  subsequent  test 
array,  or  validation  experiment,  is  performed  to  assess  the  accuracy 
of  the  model  proposed.  Both  oxidized  and  cleaned  foils  were  sub¬ 
jected  to  UC.  Oxidized  foils  are  not  expected  to  bond  while  cleaned 
foils  are.  This  is  due  to  the  thick  oxide  layer  present.  The  bonding 
process  relies  on  the  breakup  of  the  oxide  layer  to  produce  clean 
metal  surfaces  that  provide  paths  of  diffusion.  In  Janaki  Ram  et  al.’s 
(2007)  work  of  UC  of  multi-material  systems  it  is  discussed  that 
if  the  oxide  is  not  broken  apart  and  properly  dispersed  voids  will 
remain  since  the  UC  bonding  process  is  relies  on  metal-metal  con¬ 
tact.  Thus  welds  using  oxidized  foils  will  be  representative  of  pure 
slip,  friction  only  condition.  Inconsistencies  in  friction  coefficient 
trends  between  oxidized  and  cleaned  foils  can,  in  part,  be  attributed 
to  potential  deformational  heating  effects,  which  are  expected  to 
be  small.  All  foils  are  consolidated  using  a  foil-foil  setup  as  shown 
in  Fig.  8.  In  this  arrangement,  both  the  horn  and  anvil  contacting 
surfaces  are  knurled  to  prevent  foil  slip  at  the  tool-foil  interfaces. 

For  MMC  tape  welds,  a  tape-tape  geometry  was  not  used.  Stack¬ 
ing  two  unprocessed  tapes  on  top  of  each  other  amplifies  the  tape’s 
geometric  non-uniformities  resulting  in  non-uniform  tempera¬ 
tures.  Fig.  9  illustrates  the  difference  in  temperature  distributions 
across  the  sample’s  width  of  a  tape-tape  and  tape-substrate 
weld.  The  samples  compared  have  similar  (4%  difference)  aver¬ 
age  temperatures,  but  substantially  different  process  settings.  Weld 
settings  for  force,  speed,  and  amplitude  were:  1451  N,  51.3  mm/s, 
and  28.4  p,m  for  the  tape-tape  weld  and  2461  N,  19.3  mm/s,  and 
25.8  p,m  for  the  tape-substrate  weld,  respectively.  While  sufficient 
energy  can  be  delivered  to  permit  gap  closure  of  a  tape-tape  sam¬ 
ple;  this  resulted  in  tool  sticking  at  the  anvil-tape  interface.  By 
welding  to  a  substrate  (Fig.  8(b)),  energy  levels  can  be  increased  and 
non-uniformities  are  greatly  reduced.  Note  that  for  tape-substrate 
welds  a  knurled  anvil  is  not  used  to  prevent  slip,  instead  the  sub¬ 
strate  is  bolted  in  place  (as  shown  in  Fig.  6). 
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Fig.  7.  Optical  microscope  images  of  a  through  thickness  cross-section  of  MetPreg.  Images  show  left  center  and  right  views. 
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Fig.  8.  Welding  geometries:  (a)  foil-foil  with  knurled  anvil;  (b)  tape-substrate  with  bolt  down  substrate. 


3.3.1.  Contact  length,  lc 

The  contact  length  between  materials  to  be  consolidated  was 
determined  by  performing  a  series  of  spot  welds  over  a  range 
of  pressures.  Spot  welds  were  performed  by  welding  a  sample 
with  a  feed  rate  of  0  mm/s  and  amplitude  of  24  p,m  for  a  duration 
of  approximately  0.2  s.  Upon  separation  of  consolidated  samples 
the  length  over  which  friction  abraded  the  contacting  surfaces 
is  measured  and  defined  to  be  the  contact  length.  The  abraded 
area  was  slightly  larger  than  the  previously  bonded  areas.  This 
prevented  any  potential  bias  that  could  have  occurred  should  tear¬ 
ing  of  the  bonded  region  take  place.  Neither  foil-foil  welds  nor 
tape-substrate  bonds  showed  significant  dependence  on  pressure 
(see  Fig.  10).  That  is  to  say,  the  apparent  contact  area  did  not  depend 
on  pressure;  however,  the  real  contact  area  will  be  influenced  by 
the  applied  load  since  the  sample’s  surface  will  deform  until  enough 
area  is  in  contact  to  balance  the  applied  load.  The  variation  in  real 
contact  area  will  alter  the  contact  geometry  which,  according  to 
Blau's  (2009)  well  documented  book  on  friction  concepts  and  appli¬ 
cations,  is  known  to  affect  the  friction  coefficient,  which  may  be 
one  of  the  reasons  force  can  influence  pt.  For  this  model,  apparent 
area  is  used,  so  variation  in  real  contact  area  is  not  quantified.  For 
foils  and  tapes  the  average  contact  lengths  were  2.55  and  7.22  mm, 
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Fig.  9.  Comparison  of  temperature  uniformity  across  the  width  of  a  tape-tape  weld 
and  a  tape-substrate  weld. 


respectively.  The  longer  contact  length  between  the  tape  and  the 
substrate  can  be  attributed  to  several  factors  including:  narrower 
tape  width,  greater  thickness,  and  softer  base  metal. 


3.3.2.  Friction  coefficient  determination 

For  foils,  Eq.  (12)  can  be  solved  explicitly  for  pt  as  a  function  of 
temperature.  For  tapes,  pi  must  be  found  iteratively  from  inverse 
modeling.  With  a  16  sample  test  array,  the  end  result  for  both  foils 
and  tapes  is  a  series  of  16,  likely  different,  friction  coefficients;  one 
for  each  weld.  These  values  are  averaged  if  one  wants  to  report  a 
single  value,  inconstant,  or  further  analysis  is  carried  out  to  determine 
a  parameter  dependent  function,  /nRSM. 

To  determine  ptConstant  all  16  values  were  averaged.  This  results 
in  the  single  valued  constant  friction  coefficient,  inconstant,  which  is 
valid  over  the  range  of  processing  parameters  tested.  To  find  plrsm  a 
response  surface  model  (RSM)  was  generated  which  fits  a  surface  to 
the  experimentally  determined  friction  coefficients.  This  is  done  by 
fitting  a  quadratic  equation,  which  includes  all  linear  combinations 


Force  (N) 

297  874  1451  2028  2605  3182 


Fig.  1 0.  Contact  length  between  consolidated  materials  as  a  function  of  applied  load. 
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of  parameters  according  to  a  second  order  polynomial,  Eq.  (16),  to 
the  experimental  data. 

P-rsm  =  b0  +  bxk  +  bFF  +  bss  +  bxxX2  +  bFFF 2  +  bsss2  +  bXFXF 

+bXsXs  +  bFsFS  (16) 

where  the  fa’s  are  the  fitting  coefficients  and  X,  F,  and  s  are  coded 
welder  parameters;  all  of  which  are  dimensionless.  Coded  param¬ 
eters  vary  between  -1  and  1  taking  the  form: 


~  X  Xavg 

X  =  - . - - 

Ax 

(17) 

A  „  Xmax  -  xmin 

(18) 

The  credibility  of  /z,  constant  or  variable,  is  assessed  by  com¬ 
paring  the  trends  in  /z  determined  from  the  experiments  to  the 
expected  trends  from  the  literature,  Fig.  2.  Once  the  coefficient 
of  friction  is  known  a  new  test  array  was  created  to  evaluate  the 
model’s  accuracy. 

3.3.3.  Parameter  selection 

The  selection  process  for  welder  parameters  differed  for  foils 
and  tapes.  To  minimize  the  effects  of  deformation,  oxidized  foils 
were  given  the  maximum  amount  of  energy  that  did  not  result  in 
welding,  thus  allowing  for  pure  slip  throughout  the  welding  pro¬ 
cess.  The  same  parameter  array  was  employed  for  cleaned  foils  to 
allow  for  comparison.  The  oxide  layer  is  much  thinner  for  cleaned 
foils,  so  bonding  did  occur  during  consolidation.  Therefore,  cleaned 
foils  may  not  slip  throughout  the  entire  weld  duration.  The  modi¬ 
fications  in  surface  characteristics  will  alter  the  frictional  behavior. 
Blau's  (2009)  book  reports  that  surface  oxide  presence  has  a  lubri¬ 
cating  effect,  reducing  //.  Therefore,  increases  in  /z  can  be  attributed 
to  the  pretreatment  acid  removal  of  lubricating  surface  oxides.  Vari¬ 
ations  in  /z  may  also  result  from  bonding/no-slip  deformational 
influences.  While  the  additive  contribution  of  deformational  heat 
is  small,  the  prevention  of  slip  could  result  in  less  frictional  heat 
generation;  if  this  is  the  case  the  apparent  /z  (i.e.  the  /z  calcu¬ 
lated  with  this  approach)  will  be  smaller  than  the  real  /z  at  regions 
experiencing  slip. 

Tape  welding  parameters  were  chosen  to  achieve  good  bond¬ 
ing  and  uniform  heating.  A  good  bond  is  one  in  which  the  tape, 
or  foil,  breaks  in  tension  or  flexure  prior  to  fully  debonding  from 
its  mating  surface.  This  test  is  performed  qualitatively.  To  date 
there  is  not  an  established  testing  method  that  can  fully  character¬ 
ize  the  bond  strength  between  UC  materials;  though,  an  extensive 
body  of  work  on  mechanical  testing  of  ultrasonically  consolidated 
thin  foils  has  been  performed  at  Loughborough  University.  In  2003 
Kong  et  al.  began  their  method  development  of  mechanical  test¬ 
ing  of  UC  Al  6061  welds.  In  2004(a)  Kong  et  al.  continued  their 
work  with  Al  3003;  methods  such  as  lap-shear  and  peel  test  have 
been  attempted;  however  typical  failure  modes  do  not  occur  uni¬ 
formly  in  the  bonded  interface.  Samples  either  break  outside  of  the 
bonded  region,  or  tear  non-uniformly  failing  to  capture  the  stable 
debonding  load.  Nevertheless,  the  peel  test  is  the  most  commonly 
used  method  in  Loughborough’s  work  and  was  recently  employed 
in  Friel  et  al.’s  (2010)  paper  to  investigate  the  effect  of  surface 
topography  for  UC  of  Al. 

Uniform  heating  is  characterized  by  removal  of  the  thermal  edge 
peaks,  or  double  peaks,  shown  in  Fig.  9.  As  load  levels  increase, 
pressure  in  particular,  greater  levels  of  uniformity  are  achieved. 
Due  to  the  thermal  mass  of  the  horn  and  substrate,  edge  effects 
cannot  be  completely  eliminated. 

The  corresponding  parameter  ranges  found  to  fit  the  desired 
welding  criterion  specified  in  this  section  are  summarized  in 
Table  3.  Note  that  tape-substrate  welds  used  higher  levels  of  all 


Table  3 

Experimental  parameter  ranges  used  during  UC. 


Foil-foil 

Tape-substrate 

Amplitude  (pan) 

9.4-18.4 

24-34 

Force  (N) 

874-1739 

1739-2605 

Speed  (mm/s) 

87-123 

10-56 

Table  4 

Initial  LI  6  Taguchi  array. 

Sample  # 

Amplitude  (fjim) 

Force  (N) 

Speed  (mm/s) 

Foils 

Tape 

Foils 

Tape 

Foils 

Tape 

1 

9.4 

34.0 

1162 

2028 

99.2 

25.4 

2 

12.4 

24.0 

874 

1739 

99.2 

56.2 

3 

18.4 

27.3 

1739 

2028 

87.2 

56.2 

4 

12.4 

27.3 

1451 

1739 

123 

40.8 

5 

9.4 

34.0 

874 

1739 

87.2 

10.1 

6 

15.4 

24.0 

874 

2317 

111 

25.4 

7 

9.4 

27.3 

1451 

2605 

111 

25.4 

8 

15.4 

30.7 

1451 

2028 

87.2 

10.1 

9 

15.4 

24.0 

1162 

2605 

123 

10.1 

10 

18.4 

34.0 

1162 

2605 

111 

56.2 

11 

12.4 

30.7 

1739 

1739 

111 

25.4 

12 

18.4 

24.0 

1451 

2028 

99.2 

40.8 

13 

12.4 

27.3 

1162 

2317 

87.2 

10.1 

14 

9.4 

34.0 

1739 

2317 

123 

40.8 

15 

18.4 

30.7 

874 

2317 

123 

56.2 

16 

15.4 

30.7 

1739 

2605 

99.2 

40.8 

parameters  than  foil-foil  welds.  This  was  necessary  to  achieve  uni¬ 
form  heating. 

3.3.4.  Parameter  implementation 

A  design  of  experiments  approach  was  taken  to  generate  the 
UC  test  array.  Design  of  experiments  is  a  method  often  used  to 
optimize  or  predict  a  system  response,  in  this  case  temperature. 
A  four  level  LI  6  orthogonal  Taguchi  array  was  first  employed  to 
determine  the  friction  coefficient  by  means  of  the  processing  tem¬ 
perature.  A  subsequent  3-by-3  test  array  was  then  used  to  assess  the 
model’s  temperature  prediction  accuracy  and  limitations.  The  LI  6 
and  3-by-3  arrays  used  for  the  foil-foil  tests  and  the  tape-substrate 
test  were  run  according  to  the  sequences  shown  in  Tables  4  and  5, 
respectively. 

4.  Results  and  discussion 

The  thermal  results  from  the  initial  test  array  and  the  result¬ 
ing  calculated  friction  coefficients  are  presented  in  this  section.  An 
analysis  of  the  experimentally  determined  friction  coefficients  is 
detailed  by  comparing  calculated  and  literature  values.  A  subse¬ 
quent  validation  test  using  the  determined  friction  coefficients  is 
executed  to  indicate  the  model’s  accuracy. 


Table  5 

Validation  (3-by-3)  array. 


Sample  # 

Amplitude  (fim) 

Force  (N) 

Speed  (mm/s) 

Foils 

Tape 

Foils 

Tape 

Foils 

Tape 

1 

16.9 

32.5 

1595 

2172 

111 

50.0 

2 

16.9 

25.8 

1018 

1884 

99.2 

50.0 

3 

13.9 

25.8 

1018 

2461 

111 

19.3 

4 

16.9 

29.1 

1307 

2461 

87.2 

34.7 

5 

13.9 

32.5 

1307 

1884 

87.2 

19.3 

6 

10.9 

32.5 

1018 

2461 

123 

50.0 

7 

13.9 

29.1 

1595 

1884 

123 

34.7 

8 

10.9 

29.1 

1307 

2172 

99.2 

34.7 

9 

10.9 

25.8 

1595 

2172 

87.2 

19.3 
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Fig.  11.  L16  Results  -  Temperature  (top)  and  friction  coefficients  (bottom). 


4.1.  LI  6  array  -  temperature  and  friction  coefficient  results 

The  temperature  results  and  calculated  friction  coefficients  from 
the  initial  LI  6  array  are  shown  in  Fig.  11.  Sample  numbers  from 
Table  4  correspond  to  the  test  numbers  in  Fig.  11.  Tape-substrate 
temperatures  were  obtained  by  averaging  all  temperature  contours 
in  the  length  and  width  direction.  Foil  temperatures  were  obtained 
by  averaging  the  maximum  temperature  recorded  for  each  contour. 
Maximum  temperature  was  used  for  foils  because  it  is  less  compu¬ 
tationally  intensive  to  obtain.  Since  the  Biot  number  of  the  foils  is 
low  (i.e.  temperature  is  spatially  invariant),  maximum  and  average 
temperatures  are  comparable.  As  expected  there  is  a  large  amount 
of  variability  among  the  computed  friction  coefficients.  A  summary 
of  friction  values  obtained  is  shown  in  Table  6.  The  range  of  results 
determined  in  this  study  fall  within  the  same  range  as  reported  in 
the  literature  (see  Fig.  2). 

The  friction  coefficient  for  tape-substrate  welding  is  substan¬ 
tially  lower  than  that  of  the  foil-foil.  This  is  a  result  of  differing 
loading  conditions  and  material  properties.  For  foil-foil  welds, 
loads  were  chosen  specifically  to  avoid  plastic  deformation.  How¬ 
ever,  this  was  not  the  case  for  the  MetPreg  tape-substrate  welds. 
For  tape  welds  the  normal  load  alone  was  in  excess  of  the 
matrix  material’s  room  temperature  yield  strength.  Additionally, 


Table  6 

LI  6  friction  coefficients. 


Average 

Minimum 

Maximum 

Coefficient  of 
variation 

Oxidized  foils 

0.399 

0.233 

0.644 

32% 

Cleaned  foils 

0.539 

0.301 

0.843 

29% 

Tape-substrate 

0.145 

0.076 

0.211 

24% 

changes  to  the  tapes  cross-sectional  area  post-welding  also  sup¬ 
port  the  presence  of  bulk  plastic  deformation  during  consolidation. 
Bhushan’s  (1999)  comprehensive  textbook  on  tribology  presents  a 
preliminary  derivation  for  plastic  contact  of  ductile  metals  which 
suggests  /t<1/5;  this  is  much  lower  than  typical  sliding  fric¬ 
tion  coefficients.  This  explains  the  decreased  friction  coefficient 
observed  during  UC  of  tapes. 

Friction  coefficients  also  differed  when  comparing  oxidized  and 
cleaned  foils.  Oxidized  foils  have  a  larger  amount  of  aluminum 
oxide,  AI2O3,  present  on  the  surface  as  compared  to  cleaned  foils. 
Bhushan’s  (1999)  textbook  reports  self-mated  friction  coefficient 
values  for  Al  are  higher  than  that  of  AI2O3.  In  addition  to  changing 
the  surface  oxidation  levels,  cleaning  can  remove  trace  lubricants 
and  other  contaminants.  Blau's  (2009)  book  describes  that  the  more 
thorough  the  cleaning  process,  the  higher  the  friction  coefficient 
can  become. 

Another  interesting  observation  is  the  amount  of  variation  dur¬ 
ing  testing  for  each  material  pairing.  The  coefficient  of  variation,  CV, 
of  measured  temperatures,  T-CV,  was  just  over  30%  for  all  mate¬ 
rial  pairings.  For  foils,  oxidized  and  cleaned,  this  correlated  to  a 
p,-CV  of  roughly  30%,  while  tape-substrate  welds  had  a  p,-CV  of 
only  20%.  Thus,  the  friction  behavior  of  an  ultrasonically  consoli¬ 
dated  tape-substrate  interface  is  not  as  easily  influenced  by  welder 
parameters  as  a  foil-foil  interface.  This  could  be  a  function  of  the 
large  difference  in  time  scales.  Weld  times  were  on  average  more 
than  1 0  times  longer  for  tapes  than  foils.  Thus  it  is  more  likely  for  the 
friction  coefficient  to  have  fully  stabilized  (Fig.  2,  ft  vs.  N).  Fiber  rein¬ 
forcement  may  also  play  a  role  by  affecting  the  abrasion  resistance 
of  the  contacting  surfaces.  By  limiting  wear,  the  surface  properties 
could  be  less  susceptible  to  external  forces  (i.e.  welder  parameter 
settings).  An  investigation  into  the  wear  mechanics  of  UC  materials 
may  provide  more  insight  into  this  phenomenon. 
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Table  7 

P-rsm  -  response  surface  model  coefficients. 


bo 

bx 

bp 

bs 

bxx 

bpp 

bss 

bxF 

bxs 

bFs 

Oxidized  foils 

0.3885 

0.0952 

-0.1547 

0.0015 

-0.0173 

0.0561 

-0.0201 

-0.0326 

-0.0383 

0.0129 

Cleaned  foils 

0.5430 

0.0897 

-0.1771 

0.0221 

-0.0470 

0.0464 

-0.0067 

-0.0082 

-0.0094 

-0.0122 

Tape-substrate 

0.1462 

0.0238 

-0.0269 

0.0087 

-0.0077 

0.0067 

-0.0018 

-0.0040 

-0.0119 

0.0018 

4.3.1.  Variable  friction  coefficient  -  RSM 

Minitab  (2003)  was  used  to  perform  the  response  surface  model 
calculation  as  described  in  the  experimental  section.  The  RSM  coef¬ 
ficients  are  shown  in  Table  7.  The  values  from  Table  7  can  then  be 
paired  with  Eq.  (16)  to  plot  the  trends  of  the  friction  coefficient  for 
each  material  pairing  over  a  range  of  amplitudes,  forces,  and  speeds 
(Figs.  12-14). 


4.2.  Friction  coefficient  validity 

There  is  good  correlation  in  the  coefficient  of  friction  trends 
when  comparing  the  experimental  results  here  (Figs.  12-14)  to 
fretting  fatigue  experiments  done  by  Naidu  and  Raman  (2005), 
Fig.  2.  In  Figs.  12-14,  L,  M,  and  H  correspond  to  Low,  Medium, 
and  High  welder  parameter  levels,  respectively.  Low  is  the  low¬ 
est  setting  tested  experimentally,  Medium  is  the  average,  and  high 
is  the  highest.  From  Figs.  12-14,  increase  in  oscillation  amplitude 
causes  an  increase  in  /x,  while  normal  force  has  a  diminishing  effect. 
This  directly  correlates  to  the  literature  findings.  On  the  other  hand 
speed,  or  number  of  cycles,  N,  has  little  effect  on  /x.  Under  isother¬ 
mal  conditions,  speed  should  have  no  effect  on  /x,  however,  UC  is 
not  an  isothermal  process,  and  thus  lower  speeds  (higher  N)  cor¬ 
respond  to  higher  temperatures.  This  introduces  the  potential  for 
error  since  //  does  depend  on  temperature  as  illustrated  by  Zhang 
and  Li’s  (2007)  experimental  results  previously  discussed,  Fig.  2.  As 
a  result  the  friction  coefficient  for  both  foils  and  tapes  exhibit  some 
dependence  on  speed  which  may  be  attributable  to  this  effect.  Thus, 
assuming  /x  does  not  depend  on  temperature  introduces  a  small 
amount  of  error. 


Another  potential  source  of  error  is  the  presence  of  defor- 
mational  heating.  This  model  assumes  frictional  heating  only 
throughout  the  entire  welding  process.  Thus,  any  deformational 
heating  gets  intrinsically  lumped  in  to  the  friction  coefficient.  How¬ 
ever,  contributions  from  deformation  are  small  (Section  1.1)  and 
the  /x  trends  across  all  material  pairings  are  very  similar  suggesting 
the  driving  force  behind  the  thermal  development  is  also  similar. 
When  welding  oxidized  foils  it  was  made  certain  that  machine 
parameters  were  chosen  so  as  to  avoid  effects  from  the  defor¬ 
mational  regime,  in  particular  bonding.  Thus,  oxidized  foils  were 
welded  under  pure  slip  conditions,  therefore  friction  dominant. 
Since  trends  in  the  friction  coefficient  from  the  friction  dominated 
oxidized  foil  welds  are  similar  for  the  other  material  pairings,  this 
supports  that  friction  is  also  the  dominant  heating  mechanism  for 
cleaned  foils  and  tape-substrate  welds. 

4.3.  Validation 

With  both  /^constant  and  /xrsm  characterized  by  minimizing  the 
error  between  the  predicted  and  experimental  temperatures,  a 
3-by-3  array  of  the  three  distinct  values  of  the  three  processing 
parameters  (F,  s  and  7.),  as  listed  in  Table  5,  was  executed  experi¬ 
mentally  and  the  values  for  the  temperature  field  were  recorded. 
The  model  was  used  with  the  characterized  value  of  both  inconstant 
and  iirsm  to  predict  the  temperatures  and  compare  them  to  the 
experimental  values.  Predicted  and  experimental  temperatures  for 
each  material  pairing  are  shown  in  Fig.  15.  Good  agreement  was 
achieved.  A  constant  friction  coefficient  results  in  an  error  of  1 6%  on 
average.  While  a  parameter  dependent  friction  coefficient  reduces 
this  error  to  7%.  The  fitting  error  of  the  RSM  is  of  the  same  order; 
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Fig.  12.  Friction  coefficient  trends  -  oxidized  foils  -  L,  M,  and  H  correspond  to  Low,  Medium,  and  High  welder  parameter  levels,  respectively.  Low  is  the  lowest  setting  tested 
experimentally.  Medium  is  the  average,  and  High  is  the  highest. 
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Fig.  13.  Friction  coefficient  trends  -  cleaned  foils  -  L,  M,  and  H  correspond  to  Low,  Medium,  and  High  welder  parameter  levels,  respectively.  Low  is  the  lowest  setting  tested 
experimentally.  Medium  is  the  average,  and  High  is  the  highest. 


therefore,  7%  error  is  very  good  agreement.  The  average,  maximum 
and  minimum  error  percentages  are  summarized  in  Table  8.  While 
fi constant  did  not  result  in  a  high  average  error,  it  did  result  in  very 
high  maximum  error.  The  highest  errors  occurred  when  force  and 
amplitude  were  out  of  phase  with  each  other  (i.e.  one  high,  the 


other  low).  As  with  the  LI  6  array  used  in  Section  4.1,  temperature 
is  reduced  to  an  average  value  for  each  test.  This  is  the  case  for  both 
measured  and  predicted  values.  For  foil  predictions  this  is  the  only 
option  since  a  lumped  parameter  model  was  employed.  For  tape 
welds  temperature  contours  can  be  plotted  across  the  foil  width, 
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Fig.  14.  Friction  coefficient  trends  -  tape  on  substrate  -  L,  M,  and  H  correspond  to  Low,  Medium,  and  High  welder  parameter  levels,  respectively.  Low  is  the  lowest  setting 
tested  experimentally,  Medium  is  the  average,  and  High  is  the  highest. 
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Fig.  15.  Model  validation  array  (3-by-3  array):  top,  oxidized  foil-foil;  middle,  cleaned  foil-foil;  bottom,  tape-substrate.  TIR  is  the  experimentally  measured  average  temper¬ 
ature  (average  maximum  for  foils  and  average  of  all  data  for  tapes)  with  error  bars  indicating  the  size  of  the  standard  deviation,  plus  and  minus.  “T( ficonstant )"  and  'T(/zrsm)” 
are  the  analytic  (foils)  or  FEA  (tape)  temperature  solutions  using  a  constant  and  parameter  dependent  (RSM)  friction  coefficient,  respectively. 


Table  8 

Model  accuracy  within  the  range  of  parameters  in  which  the  friction  coefficient  is 
characterized.  Error  between  model  and  experiments  reported  in  percentages. 


Oxidized  foils 

Cleaned  foils 

MetPreg  L16-T  (°C) 

T(  /X  constant) 

T(I^rsm) 

Tift  constant 

)  T(^rsm) 

T(  ^constant) 

T(/Xrsm) 

Avg 

15% 

9% 

16% 

5% 

17% 

7% 

Min 

2% 

4% 

3% 

1% 

5% 

2% 

Max 

36% 

21% 

46% 

8% 

52% 

16% 

and  when  plotted  do  follow  a  parabolic  distribution  similar  to  that 
of  the  measured  data  in  Fig.  6.  To  facilitate  quantitative  comparisons 
the  average  temperature  was  used. 

5.  Conclusions 

The  processing  temperature  during  ultrasonic  consolidation  of 
thin  foils  and  prepreg  metal  matrix  composite  tapes  can  be  pre¬ 
dicted  using  a  friction  work  heat  generation  model.  Many  factors 
can  influence  the  magnitude  of  the  friction  coefficient.  Surface  con¬ 
dition  is  one  such  factor,  oxidized  vs.  cleaned.  Cleaned  surfaces 
provide  more  sliding  resistance  than  oxidized  surfaces.  Stress  state 
also  affects  the  amount  of  friction  present;  plastic  deformation  sup¬ 
presses  the  friction  coefficient.  The  friction  coefficient  can  be  taken 


to  be  either  constant,  or  parameter  dependent.  However,  literature 
and  our  experiments  show  significant  dependence  of  the  coefficient 
of  friction  on  the  machine  parameters  during  ultrasonic  consolida¬ 
tion.  Model  accuracy,  in  particular  for  unbalanced  loads  (i.e.  high 
force,  low  amplitude),  is  much  improved  when  using  a  parameter 
dependent  friction  coefficient.  The  model  developed  here  can  now 
be  applied  to  other  aspects  of  the  ultrasonic  consolidation  process 
such  as  characterizing  the  thermo-mechanical  response  or  bond 
strength  and  development. 
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